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Abstract 

We present preliminary measurements of branching fractions of the b — ► sss penguin-dominated 
decays — > (jyK^ and B° — > cj)K° in a sample of approximately 89 million BB pairs collected by 
the BABAR detector at the PEP-II asymmetric-energy i?-meson Factory at SLAC. We determine 
B(B ± -► = (10.0lg;l (stat.) ±0.5(syst.)) x 10" 6 and B(B° -> = (7.6±J;| (stat.) ± 

0.5(syst.)) x 10~ 6 . Additionally, we measure the charge asymmetry AcpiB^ —>■ (pK^) = 0.039 ± 
0.086 (stat.) ± 0.011 (syst.) and set an upper limit on the CKM- and color-suppressed decay B^ — > 
<jm ± , B{B ± -> <jm±) < 0.41 x lO" 6 (90% CL). 
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1 Introduction 

5-meson decays with a in the final state present a special interest because they are dominated by 
the b — > s(d)ss gluonic penguins, possibly with a significant contribution from electroweak penguins 
(Figs. Ei, QJ)), while other Standard Model (SM) contributions are strongly suppressed pQ. Since 
contributions of diagrams with a c or a u quark in the loop are small, all dominant SM decay 
amplitudes have the same weak phase, leading to a very small (~ 1%) predicted value of direct 
CP asymmetry Acp in — ► (pK^. However, many models of new physics introduce new heavy 
particles, with new couplings, that would contribute to these decays, potentially making Acp quite 
large The amounts of CP and flavor violation observed in these decays can therefore be used 
to constrain the parameters of models of new physics. 




Figure 1: Examples of penguin diagrams for (a) B — > <pK, (b) B — > cj>K and B — > </>7r; (c) the rescattering 
diagram for B — > <f>K and B — > </>7r. The unlabled spin-1 propagator in diagram (b) should be interpreted 
either as a hard gluon (with two or more additional soft gluons that are not shown), or as a 7 or a Z . The 
tree part of diagram (c) could lead directly to a B — * <ph final state via the small, poorly understood uu 
component of the <f> resonance 



Recent preliminary results from BABAR and Belle on the time-dependent CP asymmetry in the 
decay B° — > tpKg have raised questions about the rescattering contribution to the B° — ► (j)K° 
decay amplitude (Fig.^i). While this cannot be computed a priori, simply from the weak couplings 
it will be larger in B — > (j)7r than in B — > (j)K by a factor of roughly cot(#c) ra 4.4, where 9c is the 
Cabibbo angle jSJ. By searching for B^ — * (j)^, our analysis can constrain the magnitude of the 
rescattering contribution to B° — > 0AT°. 

Additional reasons to be interested in a detailed study of the b — > s(d)ss processes include their 
sensitivity to QCD dynamics [HI El and to the poorly measured Cabibbo-Kobayashi-Maskawa 
matrix element Vt s . 

The decays -B 1 * 1 — ► t/)^ 1 * 1 and i3° — > have previously been observed by CLEO |H1 , BABAR [0] , 
and Belle jlUj . The significantly increased size of the BABAR data set and an improved analysis 
technique allow us to achieve a substantial reduction of both the statistical and the systematic 
errors on the branching fractions of the two decays. The analysis is based on a multivariate 
maximumdikelihood fit; the yields for the decay modes B^ — > (pK^ 1 and B^ — > 07T 111 are obtained 
simultaneously. A blind analysis technique is used to avoid the potential for an experimenter- 
induced bias: the signal region is hidden until all significant details of the analysis are finalized. 
The determination of systematic errors is completed subsequently. 

2 The BaBar Detector and Data Set 

The data were collected with the BABAR detector [TT] in 1999-2002 at the PEP-II asymmetric- 
energy e + e~ collider located at the Stanford Linear Accelerator Center. An integrated lu- 
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minosity of about 82 fb was recorded at the peak of the T(4S) resonance, corresponding to 
88.9 ± 1.0 million BB pairs. 

The asymmetric beam configuration provides a boost to the T{AS) in the laboratory frame 
(/3-y k, 0.56), increasing the maximum momentum of the I?-meson decay products to 4.4GeV/c. 
Charged particles are detected and their momenta measured by a combination of a silicon vertex 
tracker (SVT), consisting of five double-sided layers, and a 40-layer central drift chamber (DCH), 
both operating in a 1.5 T solenoidal magnetic field. The tracking system covers 92% of the solid 
angle in the center-of-mass (CM) frame. The track-finding efficiency is, on average, (98 ± 1)% for 
momenta above 0.2GeV/c and polar angles greater than 0.5 rad. Photons are detected by a Csl 
electromagnetic calorimeter (EMC), which provides excellent angular and energy resolution with 
high efficiency for energies above 20 MeV. 

Charged-particle identification is provided by the average energy loss (dE/dx) in the two track- 
ing devices and by the novel internally reflecting ring-imaging Cherenkov detector (DIRC) covering 
the central region. A tt/K separation of better than 4cr is achieved for tracks with momenta below 
3GeV/c, decreasing to 2.5cr for the highest momenta arising from U-meson decays. Electrons are 
identified with the use of the tracking system and the EMC. 

3 Event Selection 

Hadronic events are selected on the basis of track multiplicity and event topology. S-meson can- 
didates are fully reconstructed from their charged decay products: <fi — ► K + K~ and K° — > K® — ► 
7r + 7r~. Charged tracks that are B or cf) daughters are required to originate from the interaction 
point (within 10 cm along the beam direction and 1.5 cm in the transverse plane), have at least 12 
DCH hits and a minimum transverse momentum of 0.1 GeV/c. Looser criteria are applied to tracks 
used to reconstruct — > it + i:~ candidates in order to allow for displaced K® decay vertices. We 
suppress e + e _ — ► t + t~ and e + e~ — > e + e~77 backgrounds by rejecting events with fewer than 5 
tracks. 

Pairs of oppositely-charged tracks that are required to originate from a common vertex are 
combined to form the <f) and K® candidates. A clean sample of K® candidates is obtained with 
requirements on the two-pion invariant mass (\M n + n - — m^o\ < 12 MeV/ c 2 ), the angle a between 
the reconstructed flight and momentum directions (cos a > 0.995) and the flight-length significance 
{I lot > 3). For (j) — ► K + K~, the invariant mass of the K + K~ pair is required to lie within the 
[0.99,1.05] GeV/ c 2 range (Fig.©. 

Tracks used to reconstruct the <p — ► K + K~ decay are distinguished from pion and proton tracks 
via a relatively loose requirement on a likelihood ratio that includes, for momenta below 0.7GeV/c, 
dE/dx information from the SVT and DCH and, for higher momenta, the Cherenkov angle and the 
number of photons as measured by the DIRC. In addition, these tracks are required to pass electron 
and proton vetoes. Determination of the flavor of the high-momentum track in B^ — ► (fch^ 
decays is provided mostly by Cherenkov-angle residuals, normalized to their uncertainties, which 
are computed for the pion and the kaon hypotheses and are used in the maximum- likelihood fit. 
During event selection, candidates are required to have Cherenkov angles consistent within ±4cx 
with either of the two hypotheses; they are also required to pass an electron veto. 

We identify 5-meson candidates kinematically using two nearly independent variables |llj : 

the beam-energy-substituted mass toes = J ((§ + Pr • Pb) 2 ) / (E^) — Pg, which is computed in the 
laboratory frame and is independent of the mass hypotheses assigned to the £>-candidate daughters, 
and the Lorentz-invariant missing energy AE = (qr'QB/ ' \fs) — \fs- Here qr and qs are four-momenta 
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Figure 2: (a) ttijck invariant mass distribution in the — > cfrh^ on-resonance sideband; (b) 
definition of the on-resonance sideband in toes an d AE 



of the T(4S) and the B candidate, s = (qr) 2 is the square of the center-of-mass energy, pr and ps 
are the three-momenta of the T(45) and the B in the laboratory frame, and Ex = qj- is the energy 
of the T(45) in the laboratory frame. For signal events, AE peaks at zero and toes peaks at the 
B mass. Our selection requires \AE\ < 0.2 GeV and toes > 5.2GeV/c 2 . Being dependent on the 
mass hypotheses assigned to the B decay products, AE provides additional momentum-dependent 
ir/K separation in the maximum-likelihood fit for B^ — > cjyh^ branching fractions. 

Detailed Monte Carlo studies demonstrate that backgrounds from other B decays are negligible. 
Backgrounds are dominated by random combinations of tracks produced in the quark-antiquark 
(qq) continuum. This background is distinguished by its jet-like structure — as opposed to the 
nearly spherical decay of the T(4S*). We have considered a variety of CM event-shape variables 
that exploit this difference. 

One such variable is | cos 9t \ , where 9t is the angle between the thrust axis of the B candidate 
and the thrust axis of the rest of the event, where the thrust axis A is defined as the unit vector 
that maximizes the thrust T = max [Y^Li \A • pl\/J2iLi \/\Pi\ 2 ^j • Since B's are non-relativistic 
in the T(4S) rest frame (/3 « 0.06), the |cos#t| distribution for true B candidates is very well 
described by a nearly flat first-degree polynomial; on the other hand, the | cos 9t\ distribution for 
B candidates found in the e + e~ — > qq continuum is sharply peaked at +1. We apply the cut 
| cos#t| < 0.9 throughout our analysis. 

Other quantities that characterize the event shape are the B polar angle 9b and the angle 9 q q of 
the -B-candidate thrust axis, both defined with respect to the beam axis, as well as the zeroth and 
the second Legendre moments of the rest of the tracks and neutrals, C n = J2Pi x L n (6i), n = 0, 
2, computed relative to the B-candidate thrust axis. For T(AS) decays into two pseudoscalar B 
mesons, the 9b distribution has a sin 2 #£ dependence, whereas the jets from continuum events 
lead to a uniform distribution in cos 9b- In 9 q g, the continuum jets give rise to a (1 + cos 2 9 gg -) 
distribution, while the thrust direction of true B decays is random. We further suppress the 
background by forming an optimized linear combination (Fisher discriminant ^Hl) of the four 
variables: | cos#b|, | cos q q\, Co and £2- 
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4 Maximum Likelihood Fit 



We use an unbinned extended maximum-likelihood (ML) fit to extract signal yields and charge 
asymmetries simultaneously. The extended likelihood for a sample of N events is 



where T>ik(xj; a) is the probability density function (PDF) for measured variables Xj of an event j 
in category i and flavor state k, and are the yields extracted from the fit. The fixed parameters 
a describe the expected distributions of measured variables in each category and flavor state. The 
PDFs are non-zero only for the correct final state flavor (k = 1 for B — > f and k = 2 for B —* /). 
In the simplest case, there are two categories, signal and background (i = 1,2). The decays with 
a charged primary daughter B^ — > cph^ (h = it or K) are fitted simultaneously with two signal 
(i = 1 for B^ — > (pK^ and i = 2 for i? 1 * 1 — ► (/ot^) and two corresponding background (i = 3,4) 



We define the event yields in each category in terms of the asymmetry A% and the total 
event yield nf. nn = riiX (l+_4.j)/2 and n«2 = JijX (1 — Ai)/2. The event yields rij and asymmetries 
Ai in each category are obtained by maximizing C. Statistical errors correspond to unit changes 
in the quantity x 2 = — 2 In (£/£ max ). The significance of a signal is measured by the square root 
of the change in x 2 when the number of signal events is constrained to zero in the likelihood fit; it 
describes the probability for the background to fluctuate to the observed event yield. 

The probability Vi{xj\a) for a given event j is the product of independent PDFs in each of 
the fit input variables Xj. These variables are AE, tweSj m KK, the Fisher discriminant J 7 , and the 
cosine of the (j) helicity angle (defined as the angle between the K + and B momenta in the (ft res t 
frame) cos Oh- In addition, in the simultaneous fit for the modes B^ — > (pK^ 1 and B^ — > fiir^ we 
include normalized residuals derived from the difference between the measured and expected DIRC 
Cherenkov angles for the charged primary daughter. Additional separation between the two final 
states is provided by AE. The AE separation depends on the momentum of the charged primary 
daughter in the laboratory frame and is about 45 MeV on average, varying from about 30 MeV for 
the highest-momentum to about 80 MeV for the lowest-momentum primary daughters available in 
our final states. If a given event has multiple combinations satisfying the selection requirements 
(which occurs in fewer than 0.2% of the events), the "best" combination is selected using a x 2 
quantity computed using all input variables with the exception, in the B^ — > (ph^ case, of the 
normalized Cherenkov-angle residuals and AE, which are used for (^^/(pK^ separation. 

The fixed parameters a defining the PDFs are extracted for signal from Monte Carlo simula- 
tion and for background distributions from the on-resonance sidebands in toes and AE (Fig. I^))- 
The MC resolutions and means are adjusted, when necessary, by comparing data and simulation in 
abundant calibration channels with kinematics and topologies similar to signal, B + — > tt + D° (D° — > 
K + tt~) and B° -> n + D- (D~ -> K°ir-). The PDFs for the Cherenkov-angle residuals are deter- 
mined from samples of D° — > K~ir + originating from D* decays. 

We employ a double Gaussian to parametrize the signal AE and mES PDFs. For the back- 
ground, a first-degree polynomial is used for AE and an empirical phase-space function ^1] is used 
for mES- The Fisher discriminant distributions both in signal and in background are parametrized 
by a Gaussian with different widths above and below the mean. The (^-resonance shape in sig- 
nal and the real-0 component of the continuum background are parametrized by the relativistic 




(1) 



categories. 
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Table 1: Summary of results. Equal production rates of B°B° and B + B~ are assumed. The total efficiency 
values include daughter branching fractions. Central values are followed by statistical and systematic errors; 
the upper limit on B(B ± — ► 4>n^) incorporates the associated systematic error. The statistical significance 
of the B ± ->■ </)7r ± signal is 0.5er. The 90% confidence interval for A C p {B ± -> 0i4T ± ) is [-0.104; +0.181] 





B ± -» 4>K± 


B ± (/>7T ± 


B° -» 


Events to fit 


14371 


2043 


Signal yield 


173 ± 15 


0.9±g;| (< 6.7 at 90% CL) 


soil 


Reconstruction eff. (%) 


39.8 


41.4 


43.5 


Total efficiency (%) 


19.6 


20.4 


7.4 


B (10~ 6 ) 


lO.Olgj ± 0.5 


< 0.41 (90% CL) 


7.6l}j±0.5 


A C p 


0.039 ±0.086 ±0.011 







spin-1 Breit— Wigner function j£] with the Blatt-Weisskopf damping factor correction convo- 
luted with a Gaussian resolution function (a = 1.0MeV/c 2 ); the combinatorial component of the 
rriKK distribution in the continuum background is parametrized with a second-degree polynomial 
(Fig. [2k). Since B — ► and B ^ cfiir are decays of a pseudoscalar particle into a vector and a 
pseudoscalar, the helicity-angle distribution for the signal is cos 2 9h] the background shape is again 
separated into contributions from combinatorial sources and from real <f> mesons, both of which are 
parametrized by second-degree polynomials with no linear terms. The Cherenkov-angle-residual 
PDFs are unit Gaussians for both the pion and kaon distributions. 

For all modes, we test the fitting procedure with background samples generated according to the 
PDFs and signal from Monte Carlo simulation, with numbers of signal and background events close 
to the expected. Signal yields are found to be unbiased. Correlations among the input variables in 
data are less than 5%. 

5 Physics Results and Systematic Uncertainties 

The results of our ML fit analyses are summarized in Tabled For the branching fractions, equal 
production rates of B°B° and B + B~ are assumed. The projections of the fit results are shown in 
Fig. |31 where we plot only a subsample of events, enhancing the signal with a requirement on the 
ratio of probabilities for each event to belong either to the signal or to the background categories. 

Systematic uncertainties in the ML fit originate from assumptions about the signal and back- 
ground distributions. We simultaneously vary all PDF parameters within their uncertainties and 
derive the associated systematic errors, which are found to be 2.0% for B(B ± — > 0K ± ), 10.9% for 
the 90% upper limit on B{B ± -» 0vr ± ), 2.8% for B{B° -» <f>K°), and 0.5% for A C p{B ± -> 0K ± ). 

The dominant systematic errors in the efficiency come from track finding (0.8% per high- 
momentum or <p — ► K + K~ track), particle identification (1% per <p — ► K + K~ track), and K® recon- 
struction efficiency (4.0%). Other minor systematic effects from event-selection criteria, daughter 
branching fractions, MC statistics, BB backgrounds and S-meson counting sum to 3.0%. Effi- 
ciency uncertainties affect the values of the branching fractions, but not their significances. The 
systematic uncertainty on Acp due to charge asymmetries in tracking and the DIRC is less than 
1.0%. 

Given the low significance of B(B il — ► 4"^^), we quote a 90% CL upper limit obtained by 
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Figure 3: Projection plots for B ± — > ^>K ± (left column) and B° — > (right column), made with a 

probability-ratio cut to emphasize the signal, for the variables (a,b) mES, (c,rfj A£", (ej) uikKi (g,h) cos9h- 
The solid (dashed) lines show the signal+background (background only) PDF projections 
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integrating the normalized likelihood distribution from zero. The limit incorporates changes by 
one standard deviation from uncertainties in PDFs and the reconstruction efficiency. 

6 Summary 

We have determined the branching fractions of the rare charmless penguin-dominated £?-meson 
decays B — > (pK^ and B° — > <ftK°, and have set a limit on the direct CP asymmetry Acp{B^ — > 
4>K^), with substantially reduced statistical and systematic errors compared to previously pub- 
lished results. The results contained in this paper are preliminary. 

The stringent upper limit on the CKM- and color-suppressed decay B^ — > fin^ provides ev- 
idence against the presence of large non-penguin or non-Standard Model contributions to the 
b — » s(d)ss decay amplitudes. 
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